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Alkaline-earth atoms in their lowest (nsnp) state are exceptionally long-lived and can be 
trapped magnetically. The nonspherical atomic structure leads to anisotropic long-range interac- 
tions between two metastable alkaline-earth atoms. The anisotropy affects the rotational motion 
of the diatomic system and couples states of different rotational quantum numbers. This paper 
develops a tensorial decomposition of the most important long-range interaction operators, and a 
systematic inclusion of molecular rotations, in the presence of an external magnetic field. This 
analysis illuminates the nature of the coupling between the various degrees-of- freedom. The con- 
sequences are illustrated by application to a system of practical interest: metastable **Sr. Using 
atomic parameters determined in a nearly- a6 initio calculation, we compute adiabatic potential 
energy curves. The anisotropic interatomic interaction, in combination with the applied magnetic 
field, is demonstrated to induce the formation of a long-range molecular potential well. This curve 
correlates to two fully polarized, low-field seeking atoms in a rotational s-wave state. The coupling 
among molecular rotational states controls the existence of the potential well, and its properties 
vary as a function of magnetic-field strength, thus allowing the scattering length in this state to be 
tuned. The scattering length of metastable ®*Sr displays a resonance at a field of 339 Gauss. 

PACS numbers: 34.20.Cf, 34.20.Mq, 31.10.-l-z 



I. INTRODUCTION 



The ability of van der Waals' equation-of-state of a real gas [1| to describe the phase transition from the gaseous 
to the liquid state is one of the manifestations of the universal importance that must be attributed to interatomic 
forces. Renewed interest in their detailed understanding has been prompted by the experimental demonstration 0,01 
of a pure quantum phase transition — Bose-Einstein condensation — in dilute gases of alkali-metal atoms. Ultracold 
conditions, which are necessary for the formation of a Bose-Einstein condensate, provide an ideal setting to observe, 
and measure precisely, weak interatomic interactions 

Alkali atoms constitute a natural choice: They possess strong electric dipole transitions from their ground state 
that lie in the visible range of the electromagnetic spectrum. This makes them ideal candidates for laser cooling 
techniques. Furthermore, as effective one-electron systems they have a magnetic dipole moment and can be trapped 
in a suitably shaped magnetic field. The experimental interest in alkali atoms triggered precision theoretical studies 
devoted to uncovering their long-range interaction properties 0, IS S 13 ■ 

Noble-gas atoms become accessible to laser cooling if they are not in their electronic ground state but in an 
excited, metastable one. Dispersion coefficients for metastable helium, 2 ^Sq and 2 ^Si, have been calculated, for 
example, by Chen |^ and by Yan and Babb |0| . The metastable states of the heavier noble gases can be written as 
(np5(n + l)s) (n = 2 for neon, and so on). The associated electron distributions are nonspherical, in contrast 
to the 2 ^5*0 and 2 '^S'l states of atomic He — and in contrast to ground-state alkali atoms. The interactions between 
such metastable nob le-g as atoms are not isotropic, and they are not pure dispersion forces. (Dispersion coefficients 
are reported in Ref. The most important new interaction is the electric quadrupole-quadrupole interaction, as 

discussed by Doery et al. Under normal conditions, anisotropic interactions between the constituents of a gas 

tend to be smeared out due to thermal averaging. Only ultracold temperatures permit their detailed study. 

Recently, alkaline-earth atoms have moved into the focus of interest. In their ground state they can be laser- 
cooled using the strong ^5*0 ^Pi transition Katori et al. exploited the extremely narrow ^5*0 — > ^Pi 
intercombination transition to achieve laser cooling of ground-state *^Sr down to a few hundred nano-Kelvin jl5j . 
Stable alkaline-earth isotopes with an even number of nucleons predominate, which therefore possess no nuclear 
magnetic moment: Their electronic spectra are free of hy perfine structure. This simplifies the interpretation of 
experimental data and facilitates comparison to theory [l^ . The analysis of the photoassociative spectrum of cold 
'^'^Ca atoms measured by Zinner et al. may serve as an illustration |l7j . 

Ground-state alkaline-earth atoms, which have a closed-shell electronic structure, do not lend themselves to pure 
magnetic trapping. Switching off the light field employed to hold the laser-cooled atoms in a magneto-optical trap, 
followed by evaporative cooling down to quantum degeneracy, is thus not an option. An interesting str ategy to 
overcome this obstacle is to work with metastable alkaline-earth atoms in their lowest (nsnp) ^P2 state [la flOj . 
which have radiative lifetimes of the order of ten minutes |2Q|| . The experimental feasibility of magnetic trapping of 
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metastable ®^Sr has already been demonstrated [2lLl23.l2^ . 

Metastable alkaline-earth atoms share with metastable noble-gas atoms (excluding He) the property of interact- 
ing through anisotropic forces. However, as was shown quantitatively by Derevianko pOj, the dominant electric 
quadrupole-quadrupole interaction is noticeably stronger among alkaline-earth atoms: Their large quadrupole mo- 
ment is a consequence of the relatively diffuse spatial distribution of the excited p-electron; the lack of spherical 
symmetry of a metastable noble-gas atom, on the other hand, is due to the hole in the compact valence shell, shielded 
by an orbiting electron of s-symmetry. Since metastable alkaline-earth atoms — in contrast to noble-gas atoms — do 
not suffer from trap losses via the phenomenon of Penning ionization ^4] , they offer the exciting prospect of quantum 
degenerate gases with pronounced anisotropic interatomic interactions. 

Derevianko and co-workers recently presented a first, simplified analysis of the ultracold collision properties of 
metastable alkaline-earth atoms j2^. They predict the existence of a long-range molecular potential well for the 
electronic state of highest Zeeman energy in a given external magnetic field. The scattering length on this potential 
energy curve can be tuned by adjusting the magnetic-field strength, thus enabling control of sign and strength of 
the effective interatomic interaction at ultracold temperatures. The numerical results in Ref. [23 are based on the 
assumption that s-wave scattering of two metastable alkaline-earth atoms in the low-field seeking state mentioned 
above is sensitive only to the rotationally invariant part of the interatomic interaction potential. However, as we shall 
show in this paper, that approximation is not justified in the case of relatively strong anisotropic coupling and does 
not lead to quantitatively accurate predictions. 

In order to provide a systematic framework for describing anisotropic interactions, we present in Secni^ detailed 
tensorial analysis of the leading interatomic interaction operators. In particular, we show that the electric quadrupole- 
quadrupole operator transforms like a spherical tensor of rank four, and we decompose the electric dipole-dipole 
dispersion operator into tensors of rank zero, two, and four. Similar tensorial ideas have been applied, for instance, 
to the anisotropic interactions between molecules [2^ |2^ |2^ and these were crucial to formally demonstrate the 
surprising existence of a vector interaction — a tensorial coupling of rank one — between a high-angular momentum 
Rydberg electron and an aniso trop ic ionic co re J29I I8(tI | . General introductions to the techniques of spherical tensor 
algebra may be found in Refs. |3l|. [s^. and |33j| . 

An important new ingredient, facilitated by our use of tensorial methods, is the inclusion of the quantum- mechanical 
rotation of an anisotropically interacting diatomic system in a magnetic field fSec. Illljl . At interatomic distances of a 
few hundred Bohr radii, electronic interaction energies and rotational energies are comparable. Anisotropic interatomic 
interaction leads to coupling between different rotational quantum states, and the classification of a diatomic eigenstate 
using just a single rotational quantum number ceases to be meaningful. 

Section Hvl is devoted to the discussion of a concrete example: metastable *^Sr. We have calculated all relevant 
long-range parameters (magnetic dipole moment, electric quadrupole moment, and dispersion coefficients), treating 
the correlation between the two valence electrons in atomic ®*Sr fully. Using these parameters, which are in generally 
good agreement with the ones presented in Ref. j^^, we obtain adiabatic potential curves by diagonalizing the 
complete diatomic Hamiltonian as a function of interatomic separation. Our calculations reproduce the phenomenon 
of a potential well in the low-field seeking state, correlating at large interatomic separations to a rotational s-wave 
state and fully polarized atoms with highest Zeeman energy. We make use of perturbation theory to shed light on 
how the interplay of the different physical operators generates the long-range molecular potential well. Our analysis 
clarifies how the magnetic-field strength affects depth and location of the potential minimum. It also shows that, at 
interatomic distances in the vicinity of and below the potential well minimum, the eigenstate acquires an appreciable 
admixture of higher rotational quantum states. Indeed, for this reason the scattering lengths we find on the basis of 
our numerical data differ substantially from those reported by Derevianko et al. |25|. 

A concluding discussion is given in Sec. \V\ We use atomic units throughout, unless otherwise noted. 

II. TENSORIAL STRUCTURE OF INTERATOMIC INTERACTION OPERATORS 
A. Leading expansion terms and Feshbach formalism 

The interaction Hamiltonian, Hi^t, of two atoms whose electron clouds do not overlap can be expanded in terms 
of inver se p owers of the interatomic distance, R, by applying techniques well known from classical electrodynamics 
[SJ, m, I33 • All operators in the resulting series can be expressed using purely atomic observables. 

Given a neutral atomic species with nonvanishing magnetic-dipole and electric-quadrupole moments, the leading 
expansion contributions up to fifth order in 1/R can be written as 



^^int = ^^dd + -ffmm + ^^dq + -ffqq, 



(1) 
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where 



■^^{xz, ■ x,^ - 3{x,^ ■ n){x,^ ■ n)} (2) 



*1,*2 

is the famihar electric dipole-dipole interaction operator. Xi^ symbohzes the position of an electron belonging to atom 
1 relative to the nucleus of that atom. Similarly, Xi^ refers to an electron in atom 2, measured in relation to the 
nucleus of atom 2. n denotes a unit vector along the diatomic axis. Using spherical multipole moment operators, 
defined as 



1i 



(fc) 



(r^j. , t^ij, , ifi^) representing spherical coordinates of vector a;,^. and Ci^rn{'&i^ , fi^ ) being related to the spherical harmonic 
Yi,m{'&ikj Vik) through the simple relation 



the electric dipole-dipole operator reads 

The basic assumption underlying this representation of Hm is that the vector n introduced in Eq. jSJ is identical to 
Ez, the Cartesian unit vector along the z-axis of a chosen reference frame. Hence, Eq. and all other equations in 
the remainder of this section, refer to a body-fixed frame. We will return to this point later when we incorporate the 
effect of molecular rotation. 

The second term on the right-hand side of Eq. describes the magnetic dipole-dipole interaction [s^ between 
atom 1 and atom 2: 

The magnetic dipole operator, 

/.W = -^B{i('=) (7) 

(/is: Bohr magneton; j'-''-': total angular momentum of atom fc; s'*^': spin of atom fc), is a spherical tensor of rank one. 
Since we restrict our treatment of magnetic effects to the dipole term, we suppress the rank index of fi^'^K In fact, at 
interatomic distances of only a few hundred Bohr radii or less, even magnetic dipole-dipole coupling is negligible. We 
include it here because it becomes the dominant interaction beyond a distance of a thousand atomic units, as we will 
show in Sec. IIVI 

The expansion term in fourth order contributing to _ffi„t is the electric dipole-quadrupole operator. 



-'?^!^'z^^il-^/3g«,^^^,«,,^^l,}, (s) 

and finally, in l/i?^, we have electric quadrupole-quadrupole coupling: 

H- - ^ , - (1) (2) , fi„(l)„(2) 

W.'iA^ + 'i'^U^}- (9) 

Electric dipole-octupole interactions, which are also proportional to are neglected. 

Our focus in this work is on heavy alkaline-earth atoms, in which relativistic effects are significant even in the 
valence shell. We will assume that the splitting between neighboring atomic energy levels corresponding to the same 
fine-structure manifold be larger than the interatomic interaction energy. Let both atoms be in the same well-defined 
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fine-structure state characterized by a total atomic angular momentum j and a collective quantum number ^ taking 
into account all other electronic degrees-of-freedom, apart from the projection quantum number m of the total atomic 
angular momentum. Without interatomic coupling, direct products of the form 

span an energetically degenerate, {2j + l)^-dimensional subspace of the full electronic Hilbert space of the diatomic 
system. This uncoupled basis of the degenerate model space formed the starting point of the approach taken by other 
authors [ill 12(1 l25l|. However, in order to make most efficient use of the powerful techniques of tensor algebra and 
also to demonstrate and exploit inherent symmetries of the problem, we prefer to work in a coupled representation: 

711 

X |j,m,0^'^ |j,M-m,0^'^ 

C{j,j, J; m, M — m, M) is a Clebsch-Gordan coefficient (39i| . mediating the transformation from the uncoupled basis to 
the coupled one with the total electronic angular momentum J running from to 2j and Af = — J, — J+1...,J— 1,J. 
The capital letter S is employed to collectively symbolize all other electronic quantum numbers of the noninteracting 
diatomic system, including for example j. 

By representing Hmt in the (2j + l)^-dimensional model space and diagonalizing the resulting matrix, approximate 
interaction energies can be obtained. Note that 

{J,M,E\Hm\J",m",E) =0 (11) 

and 

(J,M,S| iJdql J",m",S) =0, (12) 

because all atomic eigenstates comprising the model space have the same parity and, therefore, are not coupled 
by the electric dipole operator. The two atoms experience no direct electric dipole-dipole and no direct electric 
dipole-quadrupole interactions. 

However, another important ingredient is still missing: the electric dipole-dipole dispersion interaction, which is a 
consequence of coupling to electronic states outside the model space. This can be included by defining complementary 
projection operators 

P:=5^|J,M,5)(J,M,5| (13) 

J,M 

and 

Q:=l-r. (14) 

Let Hq denote the Hamiltonian of the noninteracting system, such that 

Ho\J,M,E)^Eo{E)\J,M,E). (15) 

(The eigenenergies of Hq depend only on atomic quantum numbers j and ^, not on J and M .) According to Feshbach 
|40l l4l| , any eigenenergy E of the interacting system, defined by the Schrodinger equation 

{i?0+-ffint}|*) =^|*), (16) 

can formally be computed within the model space, provided V I'i^) 0. P |^'), an element of the model space, is an 
eigenvector of the effective Hamiltonian 

iJeff := V{Ho+Hi,,t}V (17) 
+VHn,tQ [E-Q {Ho + i/int} QHi^tV; 

the associated eigenvalue is E: 



(18) 
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By making use of Eqs. iQJ, Hll() . H12() . and (|15|l . corrections to £'o(^) due to interatomic coupling can be found by 
diagonalizing 

Hi := i/eff - Eo{E)r = {i?n.m + i/qq + i/dis} V. (19) 

i/dis is the electric dipole-dipole dispersion interaction operator, which is proportional to 

Hdis HmQ [EoiE) ~ QH^QY^ QH^- (20) 

Expansion terms of higher order in 1/R (as well as magnetic dipole-dipole dispersion coupling) have been dropped in 
Eq. 

B. Coupled representation of interatomic interaction operators 

Next we address the transformation properties of Hi under simultaneous rotations of atom 1 and atom 2. Since 
excited alkaline- earth atoms are, in general, nonspherical, it is natural to ask how the interatomic interaction reflects 
this lack of isotropy. To that end we seek a compact tensorial formulation of i?mm, Hqq, and ffdis, by coupling the 
respective atomic tensor operators to diatomic ones. This approach also has the advantage that it enables a very 
systematic and elegant evaluation of matrix elements of Hi with respect to the coupled basis introduced in Eq. HlOfl. 

We first turn our attention to i?mm, Eq. (0. Each direct product /imV-m of atomic tensor operators can be 
expanded as follows: 



where 



Mko ■■= Ifi^^^ <E) fi^^A (22) 

L iK,0 

= ^C(l,l,X;m,-m,0)/iW/i^l 



is the 0-component of the irreducible tensor product of rank K of tensors /.t'^' and fi'-^\ Upon inserting Eq. H21|l in 
Eq. (O, it is easily seen that 

/fi 

ifmm = -^M2,0, (23) 

which demonstrates the magnetic dipole-dipole interaction operator is the 0-component of a second-rank tensor. A 
similar argument can be applied to show that Hqq (Eq. lO) is the 0-component of a fourth-rank tensor: 



Hqq = ^Qifi, (24) 



70. 

(25) 



4,0 



For the purpose of uncovering the tensorial structure of the dispersion interaction operator, iJdis (Eq- H20|l ). we 
draw on considerations used by Fano and Macek in their study of the angular distribution and polarization of the light 
emitted by atoms excited in collision processes . Greene and Zare exploited those ideas to describe the anisotropic 
emission of photofragments 43] and of laser- induced fluorescence ■ 

In analogy to Eq. (|23|l . the electric dipole-dipole interaction operator in Eq. H20|l is 

/r 

Hm = -^^2,0, (26) 



2,0 



(27) 



Furthermore, the operator Q [£'o(^) — QHqQ] ^ Q is a tensor of rank zero: 



Oo,o ■■= Q[EoiE) - QHoQ] ^Q. (28) 
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To prove this, we employ a more explicit representation of Oq,q: 

1 



Oo.o = E 



(29) 



X \J' ,m',e' ){j' ,m',e' 
j' ,m' 



Each operator J2m' I"^ i i ) {J ,^ | is, according to Racah's definition of irreducible tensors |43,|4a,|43, a 
scalar, i.e. a tensor of rank zero, since it is not difficult to sec that for m = +1,0, —1 



j™,E|j',m',s' )(j',m',s' 



= 0. 



(30) 



J-fi, Jo, and J_i are the spherical tensor components of the total electronic angular momentum operator J. Oo,0i 
being a sum of scalars, is therefore also a scalar. The dispersion interaction operator can now be written as 



Hdis — 



_6_ 



i?6 



D2,oOo,qD2,0 



J2 C(2,2,X;0,0,0)//^,, 



K=0 



1^ ^^/0,0-^^/2,0 + ^1/4,0 , 



the tensor operators Ik,o being given by 



[[D2 (E> Oo] 



D 



K = 0,2,4. 



(31) 



(32) 



In this way, a decomposition of the effective interaction operator Hi into a scalar, a second-rank tensor, and a 
fourth-rank tensor has been achieved: 



Hi 

To,o 

T2,0 
/4,0 



- P{To,o + T2.o + Ti^o}V, 
6 [2^ V6,^ 



(33) 



(The projection operator P is a scalar, in view of Eq. H30|) .) The physical significance of this result lies in the distinct 
transformation properties of the three different tensor types under rotations generated by J. These are rotations of 
the two atoms about their respective centers. The scalar term, Tq.o, is invariant under such rotations: It describes 
a purely isotropic interaction. The anisotropy of the long-range interaction between two metastable alkaline-earth 
atoms is a consequence of the presence of the two higher-rank tensors. r2,o transforms like the spherical harmonic 
^2,0 (a "d-orbital"), T^^ like 14,0 (a "g-orbital" ) . Of course, rotations about the z-axis — the interatomic axis — leave 
all tensorial terms of Hi invariant. 



Reduced matrix elements 



We next focus on the representation of Hi with respect to the coupled basis of the model space, Eq. pO|l . By 
putting to use the celebrated Wigner-Eckart theorem >4B\, the matrix elements of the operator M2.0, for instance, can 
be written as 



(J,M,S|M2,o| J ,M ,S) = 



1 

V2J-M 



C{J ,2,J;M ,0,Af)(J,S |1 M2 \\ J , S) 



(34) 
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The quantity (J, 2 || M2 \\ J ,2) is referred to as a reduced matrix element. It is independent of projection quantum 
numbers and, for that reason, is not specific to the body-fixed frame in which these expressions have been evaluated. 
Using the definitions 

[A, B, Z] := {2A + 1){2B + l)...(2Z + 1) (35) 

and 

U, C II /^^'^ II J, = {j, e II /^^'^ II J, =■■ {j, e II M II J, , (36) 

standard Wigner-Racah algebra allows a straightforward evaluation of the reduced matrix elements of the tensor M2'- 





j 








j 






j 






li 


j 






J" 








j" 





{ijj)J,E II ® M^'^l II S) = [J,2, {j ] 1 I (j,^ II II j^if . (37) 

^ [ J j" 2 J 

The 9-j symbol in this expression satisfies the equation 

„ { 2 j 1 1 

(38) 

This is a simple consequence of the symmetry properties of 9-j symbols under interchange of two of its rows and the 
fact that both atoms have, within the model space, the same angular momentum, j. (The two atoms would not have 
to be of the same species, however.) Hence, ( J, 2 || M2 || J , 2) = 0, if J + J is odd. In other words, states of even 
and odd J are not coupled by M2- We will show below that this is true for all operators defining Hi (Eq. (|33|l ). 
The reduced matrix element (j,^ || /i || j, ^) in Eq. (|37|l can be related to the atomic magnetic dipolc moment, 

(m> (i,J,ClMo|j,J,e) (39) 

o-,e iimII 



(j + l)(2j-|-l) 



so that 



(J, 2 II M.IU", 2)= 751 J,J"]V2 , , A il+2m±2l <^^f , (40) 

[ J j" 2 } ^ 

Equation 1)41) |l is valid only if j is different from zero. Otherwise, the reduced matrix element vanishes, because in this 
case the first and second row of the 9-j symbol in Eq. H37|l do not satisfy the triangle condition and the 9-j symbol 
therefore vanishes. 

Finding the reduced matrix elements of Q4 is equally easy. We introduce the atomic quadrupolc moment as 

(92) 2(j,j,ek2,o|j,J,0 (41) 



= V(^- + l)(SVlJ(2,+3) ^^■'^'"'^^ '"^'^^ 



and obtain for j > 1 



Because of properties of the 9-j symbol, the reduced matrix element ( J, 2 || (54 || J , 2) vanishes if j = or 1/2 or if 
(— 1)'^+'' is negative. 
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Evaluating the reduced matrix elements of the dispersion interaction tensors Ik {K = 0, 2, 4) is slightly more 
involved. We apply the well-known rules for determining reduced matrix elements of an irreducible tensor product of 
noncommuting tensor operators and arrive at 

(J,S II Ik II J", S> = E E { J" J / } (43) 

^^ (J,S|| D2 II J',S')(/,S' \\D2 II J",S) 

£;o(s)-£;o(s') 

The coupling by the interaction tensor of rank K between the two electronic states of angular momentum J and J , 
respectively, is controlled by the 6-j symbol 5^ in Eq- (|43ll . In particular, K, J, and J must satisfy the triangle 
condition \ J — K\ < J < J + K. The reduced matrix elements of D2 require a treatment similar to the one that led 
to Eq. (|37|l . Before writing down the final result for (J, S || || J 1 , however, we need to discuss the notation we 
will employ. 

Let e{j^''\ j,Q denote the transition energy associated with the transition of atom k from the state \ to 
I ) , and f {j^^\£,''^\i,i) be the corresponding dipole oscillator strength |50| : 

f{f'\e'\j,i) (44) 
3 2j + 1 

X(j,^ II qi ||/'^\e^'=M(j''\^^'^ II 91 IIJ,0- 
We use these to define an intermediate dispersion coefficient 

(-l)i+^<^'-^'^' (45) 

^>r-'>r-' /(J(^^c^^^J^c)/(J•^^^^^^^J,c) 

^ ^ h + £(j , e'~> , j: } e(j , , J, cmj^^^ , , J, ■ 

^j(fc) indicates a sum over all atomic eigenstates with angular momentum excluding ^'■-'^^ — ^ if j'*-''-' — j. The 
relation between our iJ-coefficients and the intermediate Cg-coefficients used by Derevianko et al. |25l | is 

Cr^-'^'=-|(2j + l)2B(,W,j(2),j,C)- (46) 

On physical grounds we expect only those j^'^^' to play a role that are consistent with dipole-allowed transitions. 
This fact is expressed by the triangle conditions that must be satisfied in the function 



^(J,j",X,j,/i),j(2)):=^(2j' + I)| 




(47) 



Specifically, j, and 1 must form a triangle for A not to vanish. 

Using definitions (|45|l and H47|) . the reduced matrix elements of Ik assume a rather compact form: 

4^ 

(J,S||/k|| J ,S) = -(2j + I)2[J,if,J (48) 

It is interesting to see why (J, S || Ik \\ J , S) vanishes if J + J is odd. First note that B{j^'^\ j,^) is 
symmetric under interchange of j^^-* and since S}-^^ and ^^^^ are just dummy summation variables, which may be 
renamed. (This statement relies on the premise that both atoms are identical and described by the same unperturbed 
quantum state | ■) An immediate consequence of this observation and the dipole-selection rules is the number of 
independent B-coefficients: For j > 1 there are six, corresponding to the combinations {j'^^^ = J + IjJ^^'' = j + 1}: 
{J« ^ J + hj^'^ = j}Aj}]^ = J + hj'^'^ = J-1}Aj^'^ = J,j'^'^ = j}Aj^'\= J.j'^'^ ^ J-i}. =j-l,j(2) = 

j — 1}. Second, by exploiting the symmetry properties of the 9-j symbols in Eq. 147|) . it can be shown that 

A(J,/',X,j,/2),^-(i)) ^ (-1)^+^"a(J, j",X,j,/i),j(2)). (49) 
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Hence, 

^ ^ A( J, j" , K, J, J , , J , J, C) (50) 

j(l) j(2) 
j(2) j(l) 

j(l) 

Ej(i) EjP) A{J,J" ,KJ,j'-^\j'-^y)B{j'-^\j'-^\j,^) in Eq. gHJ is therefore identical to zero if J + j" is odd. 

III. MOLECULAR ROTATIONS IN EXTERNAL MAGNETIC FIELD 
A. Invariant formulation of interatomic interaction Hamiltonian 

Suppose each of the two atoms were in a spherically symmetric electronic state (j = 0, m = 0), defined with 
respect to a quantization axis held fixed in the laboratory frame. In this case, the interatomic interaction energy 
clearly does not depend on the orientation of the diatomic axis relative to the quantization axis. If, however, the 
atoms are nonspherical (j > 0), then the anisotropic forces, as expressed through the tensorial structure of the effective 
interaction operator Hi in Eq. (|33|l . attempt to align the atomic angular momenta along the diatomic axis. The relative 
orientation of diatomic axis and laboratory-fixed quantization axis now matters. This leads to a coupling between 
electronic degrees-of-freedom and molecular rotations. In order to put this intuitive picture into more quantitative 
terms, we need a formulation of Hi that makes use only of observables specific to the laboratory frame. 

Let J7 be a rotation operator from the laboratory z-axis to the body-fixed z-axis. This operator connects each 
tensor T^fl in Eq. with the 0-componcnt of the tensor Tk measured in the laboratory frame: 

T^^f ) := U-'TKflU. (51) 

On the other hand, the operator UT^q^^U'^^ is just a linear combination of tensor components T^^\ M = 
—K, +K. The expansion coefficients are the elements of the appropriate rotation matrix: 

rA-,o=E^S(^>^>0)4^M . (52) 

M 

Here, ?? and ip are the polar and the azimuthal angle, respectively, of the diatomic axis in the laboratory frame. Since 
for integer K 

2?g(^, 0) = i-l f ' Ck^-m{^, ^) (53) 

(see, for example, Ref. j39|), the expansion coefhcients in Eq. (|52|l can also be interpreted as the components of a 
spherical tensor of rank K, defined with respect to rotations generated by L, the angular momentum operator of the 
diatomic axis. Thus, 

TKfi = V2K+1 \Ck ® Tk]o_o , (54) 

which demonstrates Tk^ is a tensor of rank zero with respect to rotations generated by the total angular momentum 
Jxot = L -\- J . Under rotations of all particles — electrons and nuclei — Hi is therefore totally invariant: 

Hi= V2irTT[CK®7'rK7']oo. (55) 

Let \L, Ml) denote an eigenstate of and L^- Eigenstates of J^^^ and Jtot.z are then generated by coupling nuclear 
and electronic angular momentum states as usual: 

|Jtot,Mtot,i, J) := ^C(L, J, Jtot;ML,Mtot-Mi,Mtot) (56) 

Ml 

X |L,Mi)|J,Mtot-Mi,S). 
We next determine the representation of Hi in this basis. Application of the Wigner-Eckart theorem gives 
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( Jtot, Mtot, j| i/i I Jtot, L , J ) = '7' ( Jtot, J II iJi II Jtot, i , J ) . (57) 

\J IJ tot + J- 

Of course, for the scalar tensor H\, Jtot and Mtot are good quantum numbers. 

L and J, however, are typically not conserved. This is most easily seen by inspection of the reduced matrix clement 
in Eq. 

{Jtot,L,J \\Hi\\Jtot,L' ,J" ) - (_i)^-t+J+L"^2Jtot + l (58) 

X V I ^, ^" ^ 

K=0,2,4 ^ ^ -^'"t 

X (i II Ck II i" )(J,S II II J",S). 
The reduced matrix element of Ck is that of a renormalized spherical harmonic (cf. Eq. Q) and is well known ^ 



(i II Ck II i" ) = VSL" + 1C(l", if, L; 0,0,0). (59) 

Hence, it vanishes if i + i is odd. Otherwise, for ii' = 2 and 4, there can be coupling between rotational states of 
different quantum numbers L and L , provided L, L , and K form a triangle. This coupling exists only if the atomic 
angular momentum j is larger than zero, because only then do (J, S || T2\\ J , S) [j > 1/2) and (J, S || 14 || J , S) 
(j ^ 1) in general not vanish. The reduced matrix elements of Tk for K — 0,2,4 follow immediately by combining 
Eqs. jnSl, l|4t)|l . (|42(l . and H48(l . and from our analysis in Sec lII ("I we already know that Hi conserves (—1)"^, but not 
the total electronic angular momentum J itself. 



B. Zeeman operator and rotation energy 

Precise experimental studies of long-range interactions can be carried out with cold, magnetic atoms trapped by 
means of suitable magnetic field configurations |5lj| . The presence of an externally applied magnetic field B breaks 
the local rotational invariance, and Jtot is no longer expected to be a good quantum number. To model this situation, 
we make the simplifying assumption that B is homogeneous: B = Be^- The complete Hamiltonian describing the 
relative motion of two interacting atoms, with reduced mass TOrcd, exposed to a magnetic field is thus 

The internal energy, £'o(S), of the two atoms (see Eq. H15|l 'l has been set to zero. 

By representing H in the coupled basis defined in Eq. (|56|) . a natural diabatic representation is obtained that 
can be used for rigorous scattering calculations. Alternatively, the sum of the matrices of the rotational kinetic 
energy operator, r'^ ' interatomic interaction operator, iJj, and the Zeeman operator, — (/x^^-* + /x^2^)o/B can 
be diagonalized as a function of R. This approach yields adiabatic potential energy curves. The dynamics on these 
curves is driven by the radial kinetic energy operator, — ^ {^W^ TiWi\' 

The calculation of the matrix elements of ftt is a trivial matter, 

II II II II L{L + 1) 
iJ^ouMtot,L,J\^-—^\J,^,,M,^,,L , J ) = ——^Sj^^^^jiJ^^^^^^ (61) 

and the treatment of Hi was the subject of the previous subsection. The Zeeman operator is discussed in the following. 
Again, we utilize Wigner-Racah algebra by noting that — (/x^^-* + (i''^^)oB is the 0-component of a rank-one tensor. It 
is independent of the orientation of the diatomic axis and therefore does not couple rotational states of different L. 

Interestingly, at least in the case of identical atoms, it also does not couple electronic states of different J. To 
understand this, consider the electronic reduced matrix elements of fi^^^ + fi^^'^ {j > 0): 

(62) 
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TABLE I: Parameters needed for a quantitative characterization of the long-range interaction between two metastable 
(5s5p) strontium atoms, (/i) (Eq. is the magnetic dipole moment, (52) (Eq. I41ll 'l the electric quadrupole mo- 

ment, and the coefficients j(2) (Eq. 1451 1 are needed to describe electric dipole-dipole dispersion coupling. The subscripts 
indicate that in the dispersion interaction process one atom makes a virtual dipole transition from j = 2 to j'^' and the other 
one from j = 2 to j^^^ . The parameters obtained in this work were calculated within a semiempirical approach. Derevianko et at 
l25l presented ab initio results; their intermediate dispersion coefficients have been converted using Eq. 1461 . The one-standard 
deviation uncertainties cited in Ref. |25| are given in parentheses. All data are in atomic units. 





This work 


Ref. [2^ 




-3.00 




{12) 


15.4 


15.6(5) 


Bi,i 


-132 


-158(16) 


B2,l 


187 


203(20) 


B2,2 


-266 


-264(26) 


Bs,i 


-343 


-415(42) 


Bz.2 


497 


555(56) 


Bi,z 


-1020 


-1290(130) 



Owing to the factor (— 1)'''+'^ -I- 1, the reduced matrix element vanishes if J -f J is odd. In addition, however, the 
6-j symbol in Eq. Ht)2|l differs from zero only if|J— 1|<J <J+1. Since, for atoms of the same species, J and J 
are integers, it can be concluded that J must equal J . The physical origin of this selection rule is the tensorial rank 
of the magnetic dipole operator. 

Using this result, it is now straightforward to write down the matrix elements of the Zeeman operator: 

( Jtot, Aftot, J| - (/x(^) + M^'^)oS I CmI,,l",j" ) = (-1)'''°'+^+''^m.„„ac/l,l"'5j,j" (63) 

xC(Jt'lt, 1, Jtot; Mtot, 0, Mtot)\/2Jtot + 1 1 / r I ^"^'^ H A^*'^ + M^'^ II J.^)B. 

* L -"tot "^tot J 

Thus we have shown that the matrix representation of the Zeeman operator is diagonal in all angular momentum 
quantum numbers except Jtot. The Clebsch-Gordan coefficient (as well as the 6-j symbol) in Eq. (|63|l imposes a 
restriction on the total angular momentum quantum numbers that can be coupled: Jtot and J^^^ may not differ by 
more than one unit. 

IV. APPLICATION TO METASTABLE STRONTIUM 
A. Atomic parameters 

In this section we present an application of the theory developed in Sees. and IIIII to collisions between cold, 
metastable ®*Sr atoms. The mass of ®*Sr [E^l is 1.60280 x 10^ a.u. The other atom-specific parameters needed are 
the magnetic dipole moment, (/i) (Eq. (|2H), the electric quadrupole moment, (52) (Eq. (|41|) 'l. and the intermediate 
dispersion coefficients B{j'-^\j^-^\j,^) =: Bj(i) j(2} (Eq. (jlSjl ). Here, j = 2, and ^ is specified in terms of the 
other quantum numbers characterizing the (5s5p) state of atomic strontium. We have performed semiempirical 
electronic-structure calculations to determine (/i), ((72), and 5^(1) ^(2). 

Atomic strontium can be regarded as a two-electron system, at least on the excitation energy scale of just a few 
electron volts relevant to the present study. This realization can be exploited to efficiently compute and successfully 
reproduce the amazing complexity of photoabsorption Rydberg spectra of alkaline-earth atoms. A detailed and general 
description of the underlying concepts, formal developments, and applications is given in Ref. | i53j . 

In a first step, we concentrate on the one-electron physics of Sr+. The single valence electron moves in the field of 
a closed-shell ionic core. The effective, spherically symmetric potential the electron experiences is assumed to be of 
the form 



Vi{r) ^^-{2 + {Z- 2) exp (-a^.ir) + a,,2r exp (-ai^sr)} " 2^ {l " exp [-[r/rif]} . (64) 
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TABLE II: Order-of-magnitude estimates of the energies, in atomic units, associated with magnetic dipole-dipole couphng 
(-Emm), electric quadrupole-quadrupole coupling (-Eqq), electric dipole-dipole dispersion interaction (iSdis), and molecular rota- 
tion (Erot) of two metastable strontium atoms separated by a distance of R Bohr radii. 



R 






Edis 


Eiot 


10 


10-^ 


10-^ 


10-^ 


10-^ 


10^ 


10-" 


10-^ 


10-3 


10-3 


10^ 


10-" 


10-" 


10-" 


10-" 


lO'' 


10-" 


10-" 


10-2^ 


10-" 



This one-electron potential is physically well motivated: At large distances from the ionic core, the electron feels the 
attraction by a practically pointlike Sr^"*" ion. As the electron comes closer, the ionic core responds to the presence 
of the electron and becomes polarized, acp = 7.5 a.u. is the dipole polarizability of Sr++ j^. Below the empirical 
cutoff radius n, there is a transition, mediated by the term proportional to a;, 2, from the exterior to the interior 
region of the ionic core. The electron interacts with an unscreened nucleus of charge Z — S8 at r much smaller than 
1/ai^i. Note that the parameters ai^i and n are assumed to be dependent on the orbital angular momentum quantum 
number I of the valence electron. They are determined by comparing the theoretical excitation spetrum computed on 
the basis of Eq. (|64(l with experimental data on Sr+ ; the values used in this work were taken from Ref . ■ 

Another important interaction operator must be included for a quantitative description of heavy alkaline-earth 
atoms: spin-orbit interaction, which we use in the form |55| 



^(so) _ s-lldVi 
' 2c2 r dr 



1-A 
2c2 



(65) 



The factor in brackets counteracts the r ^-divergence of r ^dVi/dr near the nucleus, s in Eq. (|()5|l is the spin of 

the valence electron, and c is the speed of light, v/'^"' was adapted by Ref. [H^ into the framework of eigenchannel 
R-matrix calculations and utilized to compute precise photodetachment spectra of the heavy alkali-metal anions. 

The quantum-mechanical motion of an electron exposed to the potential Vi + vj^^""^ is solved using a finite-element 
basis [H^IHsllH^ for the radial degree-of-freedom. A radial box size of 25 Bohr radii and 600 finite-element functions, 
corresponding to a quadratically spaced grid of 200 sectors inside the box and three functions in each sector, have 
proved to be sufficient. The Sr'^ model Hamiltonian h, which comprises kinetic energy and effective potential, is 
represented and diagonalized in a product basis of finite-element and spin-angular ;48i] functions up to / = 6 (i-orbitals). 
For each combination of orbital and total angular momentum quantum numbers of the single valence electron, the 18 
energetically lowest radial eigenfunctions — excluding of course the filled inner-shell states — are selected. A two-electron 
basis set is then constructed by coupling pairs of the selected eigenstates of the one-electron Hamiltonian /i, performing 
the angular momentum algebra within the jj-coupling scheme. Denoting the electron-electron Coulomb repulsion as 
l/''i2j the eigenenergies and eigenstates of the valence shell of atomic strontium are obtained by diagonalizing the 
matrix representation of h(l) + h{2) + l/ri2 in the two-electron basis set. In this way the complicated dynamics of 
the two correlated valence electrons is treated to a high degree of accuracy. 

The direct numerical evaluation of Eqs. H39|l . H41|) . and (|45|) using the two-electron eigenvectors is straightforward. 
Our results for the atomic parameters, which are converged with respect to all basis set parameters described above, 
are shown in Table together with the ab initio data presented by Derevianko et al. p5l | . The magnetic dipole 
moment of metastable ®®Sr was not explicitly discussed in Ref. I^^. In a state it seems natural to assume an 
atomic spin of 1, consistent with a magnetic dipole moment of —3. However, due to spin-orbit interaction the total 
spin is not really a good quantum number, and it had to be checked that the degree of spin symmetry violation is 
negligible. This is, in fact, the case. Our calculated electric quadrupole moment agrees well with that of Ref. [25l |: 
the intermediate dispersion coefficients also agree, to within the 10% one-standard deviation uncertainty of the data 
quoted by Derevianko and co-workers p5l |. 



B. Diatomic potential energy curves and scattering lengths 



Before we turn our attention to quantitative results on the interaction energies between metastable Sr atoms, it is 
helpful to explore the relative importance of the different interaction operators we have analyzed in Sees. ITU and IIIII 
To that end we have collected in Table UTI simple order-of-magnitude estimates as a function of interatomic separation. 
They are easily obtained by using the parameters in Table ^ and diagonalizing the individual terms the interatomic 
interaction operator. Hi (Eq. (|33|l ). consists of, represented in the coupled basis of the electronic model space, Eq. (|10ll . 
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FIG. 1: Adiabatic potential energy curves of two metastable strontium atoms ((5s5p) ^P2) in a magnetic field of 100 G, 
obtained by diagonalizing the matrix V defined in Eq. (Ifcifcill . The total angular momentum projection quantum number, Mtot, 
is +4. The marked curve correlates at large interatomic separations to two fully polarized strontium atoms (atomic projection 
quantum number m = +2) in a rotational s-wave state (L = 0). A basis describing molecular rotations in the laboratory frame 
with elements up to L = 40 has been employed in the calculation. 

Table im allows us to draw several important conclusions. Between distances of 100 and 1000 Bohr radii the electric 
quadrupole-quadrupole interaction dominates. Only as i? = 10 is approached, the electric dipole-dipole dispersion 
interaction becomes comparable. However, at these relatively small distances the interaction energies are comparable 
with the fine-structure splitting between the (5s5p) and the (5s5p) "^Pi state, which is of order 10~^ Hartree. 
In other words, at distances much lower than 100 Bohr radii the electronic model space we have chosen in Sec. Ullis 
no longer appropriate. Truncating all interaction operators of orders higher than six in 1/R (Eq. also becomes 

highly questionable under these circumstances. At i? > 100 a.u., however, the formalism developed in this paper may 
be expected to be useful for quantitative predictions. 

In the long-range limit, at R larger than 1000 Bohr radii, magnetic dipole-dipole coupling remains the only relevant 
interatomic interaction mechanism. Also shown in Table ^ are the characteristic energy quanta associated with 
the rotational motion of the diatomic system. They are comparable with dispersion energies around distances of 
R — 100, but become dominant at interatomic separations of 1000 Bohr radii and more. As is well known, at 
ultracold temperatures only rotational s-waves remain unaffected by the long-range rotational barrier and can probe 
interatomic interaction properties. In the presence of an external magnetic field, another energy scale must be taken 
into consideration: the Zeeman splitting, which is of order 10"^ Hartree for a magnetic field strength of 100 Gauss. 
Hence, typical laboratory fields are unable to distort the spin-orbit coupling pattern in ^^Sr, and our assumption of 
uncoupled fine-structure states remains valid. 

Figures ^ 121 a-^d IHI display adiabatic potential energy curves of two metastable strontium atoms exposed to a 
magnetic field of 100 G. The curves were calculated by diagonalizing, as a function of R, the matrix 

(Jtot, Aftot, L, J\ +Hi- (/x(i) + I j1, aC l", j" ) (66) 

= : (Jtot, Aftot, ^1 ^ I •/tot: Aftot, L ), 

i.e., the sum of the matrices of the rotational kinetic energy operator (Eq. H61I) ). the interatomic interaction operator 
(Eqs. l|57|l and l(55|l '). and the Zeeman operator (Eq. Exploiting the symmetry properties of V, we selected a 

basis characterized by even rotational and even electronic angular momentum quantum numbers. In our numerical 
studies we included rotational quantum numbers up to Lmax = 40, which gave converged results. 

Of immediate interest for magnetic trap experiments are low-field seeking atomic states. For j — 2 these are the 
ones with projection quantum numbers m = +1 and +2. If both atoms are fully polarized (m = +2) and in a 
rotational s-wave state (this is associated with a molecular eigenstate only at large interatomic separations), then 
the total angular momentum projection quantum number, Aftot, is +4. Aftot is conserved at all stages of a collision 
between the two atoms. The energetically highest Zeeman manifold with Aftot = +4 is depicted in Fig. ^ The 
potential energy curve correlating to the rotational s-wave is indicated in the figure. The behavior of the s-wave 
energy differs significantly from that of the curves deriving from higher rotational states, which are clearly repulsive. 
The s-wave curve is attractive at distances larger than 150 Bohr radii and becomes repulsive at smaller radii. Near 
the minimum of the resulting potential well there are several avoided crossings, which lead to inelastic losses through 
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FIG. 2: Adiabatic potential energy curves of two metastable strontium atoms ((5s5p) ^P2) in a magnetic field of 100 G. In this 
case Mtot — +3. The energy region shown corresponds to the Zeeman manifold deriving from the situation that one atom is in 
an m = +2 state and the second one in an m = +1 state. The atomic projection quantum numbers lose their validity below 
distances of a few hundred Bohr radii, where atom-atom interactions become important. 



nonadiabatic transitions to lower lying Zeeman states. Because of the steepness of the diabatically crossing channels, 
however, we do not expect severe losses and neglect them in this paper. 

If one of the colliding atoms has m = +2 and the other one m = +1 at large distances, then Mtot = +3 assuming 
L — Q. The final combination of two low-field seeking atoms — m = -1-1 for both atoms — corresponds to Mtot = +2, 
again implying an s-wave collision. The two cases, Mtot = +3 and Mtot = +2, are shown in Figs.|21and|31 respectively. 
The most dramatic difference between the cases Mtot = +4 (Fig. and Mtot = +3 (Fig. ^ is the strength of the 
attractive potential an s-wave experiences: For Mtot = +3, this is stronger by about one order of magnitude. In this 
Zeeman manifold the couplings are so strong that even a state correlating to a higher partial wave is subject to strong 
attractive forces. In Fig.|31 Aftot — +2, the potential energy curves are even more complicated and suggest interesting 
cold-collision dynamics. There is a double degeneracy in this energy region at large interatomic separations, since 
two atoms with m — +1 have the same energy as one atom with m = +2 plus a second atom with m = 0. The latter 
atom cannot be magnetically trapped. 

In the following we are going to concentrate on the long-range potential well in the energetically highest Zeeman 
manifold with Aftot = +4 (Fig.^). Molecular rotations are found to play a key role in its origin. We have calculated 
adiabatic potential energy curves for three different rotational basis set sizes and extracted those curves that correlate 
to two fully polarized strontium atoms in a rotational s-wave state. The results are shown in Fig. ^ If in addition 
to s-waves only d-waves are included — imax = 2 — the potential energy curve is purely attractive (the same holds 
true if a pure s-wave basis is used) . Only as soon as g- waves are taken into account (imax — 4) is there an effective 
repulsion below R = 150, thus leading to the emergence of a potential well minimum at large R. Nevertheless, the 
structure of the curve is still not converged. The potential energy curve obtained using L,„ax = 8 illustrates the need 
for i- and k-waves. Larger rotational basis sets make it difficult to isolate a smooth potential energy curve due to 
the appearance of pronounced avoided crossings (see Fig. P). However, the structure of the effective diabatic curve 
correlating at large distances to an s-wave is essentially the one found for Lmax = 8. 

Deeper insight into the formation of the long-range potential well and the role played by molecular rotations can 
be gained by resorting to a perturbative approach. For that purpose we take the diagonal of the matrix V defined 
in Eq. to represent the unperturbed problem; the perturbation is given by the off-diagonal elements. At large 
interatomic separations the state vector of two fully polarized metastable strontium atoms in a rotational s-wave state 
reads 



Up to first order its energy is 



|Jtot,Mtot,i, J) = |4,4,0,4) (67) 
l*o). 



(<fo|^|<i>o)=-^- 2(^)6, (68) 



where Ce.o = 5.33 x 10'^ a.u. Hence we see that the long-range attraction in the s-wave channel is entirely due to 
dispersion forces; more precisely, it derives from the scalar tensor Tq (see Eqs. (|58f) . and H59|) ). 
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In order to evaluate the second-order correction, all basis vectors must be determined that can couple directly 
to \^o). Since |$o) is an eigenvector of both the rotational kinetic energy operator and the Zeeman operator, the 
coupling must be mediated by the anisotropic interatomic interaction. Hi preserves Jtot and couples |$o) to the basis 
vectors 



l*l> 


|4,4,2,2), 


1*2) 


|4,4,2,4), 


1*3) 


|4,4,4,0), 


1*4) 


|4,4,4,2), 


1*5) 


|4,4,4,4). 



(69) 



The first two are rotational d-waves, the other three are g-waves. Their first-order energies, together with that of 
|^>o), are depicted in Fig. |31 Note that the basis vectors | Jtot, -^tot, >/) are eigenvectors of the Zeeman operator 
only if L or J vanishes. Therefore, the long-range limits of the first-order energies do not in general agree with the 
Zeeman energies of two noninteracting atoms. In view of Fig. \^ it is tempting to conclude that the strong repulsion 
that gives rise to the molecular potential well is due to the d-wave |$2), whose energy comes close to and eventually 
crosses the energy of the s-wave channel. 

However, this is not true. In Fig.Elwe have plotted the modulus squared of the coupling matrix elements ($,;| V \ ^o), 
j = 1, ...,5. Note the use of a logarithmic scale along the ordinate. The coupling of the s-wave to the g-waves {^3), 
|$4), and |<I>5) turns out to be stronger by many orders of magnitude than the coupling to the d-waves j^i) and 
1^2) ■ The reason is the tensorial structure of Hi (Eqs. and (|55|l '): The d-waves couple to |$o) through the 
relatively weak second-rank tensor T2; the g-waves, however, couple to the s-wave via the fourth-rank tensor T4, 
which is dominated by the strong electric quadrupole-quadrupole interaction. Because [$5) is closest in energy to 
l^o) (Fig. 0) and its coupling matrix element is by far the greatest (Fig. O, it is likely that it accounts for most of 
the second-order correction to the s-wave energy. Figure [7| confirms this. The difference between the full second-order 
energy and that due to [$5) alone is relatively small. What is even more important is the difference to the first-order 
energy. In first order, the s-wave experiences pure attraction by scalar dispersion forces. As a consequence of strong 
quadrupole-quadrupole coupling to g-waves associated with energetically lower Zeeman manifolds, a steep repulsive 
potential wall appears. This requires that coupling is introduced at least to second order in perturbation theory. 

The simple perturbative approach does not provide quantitative agreement with our numerical results, but it 
evidently captures the basic physics underlying the existence of the long-range potential well. Moreover, it predicts 
that the properties of the potential well depend on B, the magnetic-field strength. The energy difference between 
|$o) and the most important g-wave, {^5), is approximately independent of interatomic distance (Fig. 13)). This can 
be exploited for a compact approximate representation of the second-order energy of I'l'o), valid at magnetic fields of 
about 100 G and higher: 

E = -2{fi)B (70) 

C6,0 1 K 

i?6 (/x) B ■ 
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FIG. 4: Dependence of the potential energy curve of the "s-wave" in Fig.Qon the largest molecular rotation quantum number, 
imax, taken into account in the construction of the matrix V (Eq. 1661 '). Only if g- waves are included does a potential 
well appear. The cause of the pronounced admixture of higher partial waves is the relatively strong, anisotropic quadrupole- 
quadrupole interaction. At large R, the curves converge to 40.502 mK, corresponding to the Zeeman energy of two metastable 
strontium atoms with m — +2 in a magnetic field of 100 G. 



The quadrupole-quadrupole coupling parameter k is 1.44 x 10^ a.u. The location, with respect to R, of the potential 
well minimum described by Eq. (|70|l scales as B~^^'^, while the well depth scales as B^^^. As the magnetic-field 
strength is increased, the effect of electric quadrupole-quadrupole coupling at a given distance R diminishes. In order 
to compensate for the larger Zeeman splittings, it is necessary to go to shorter distances. Along the way the dispersion 
attraction grows stronger, until it is eventually overwhelmed by repulsion induced by electric quadrupole-quadrupole 
interaction. The potential well therefore gets deeper and its minimum is shifted to smaller interatomic distances as B 
is ramped up. 

The tunability of the molecular potential well could prove useful for controlling the properties of metastable, mag- 
netically trapped ®^Sr. The most important quantity in this context is the scattering length a fisj, which characterizes 
the effective interaction between ultracold collision partners |60j |. A positive a describes effective repulsion; negative 
a implies attraction. The strength of the interaction is governed by |o|. 

We have therefore calculated the scattering length in the long-range potential well as a function of magnetic-field 
strength, B. Our results are displayed in Fig. |S1 For each B value, we determined, in a procedure similar to the one 
that led to Fig. 01 the maximum L that gave a smooth potential energy curve compatible with the respective result 
for Lniax = 40, but which no longer exhibits pronounced avoided crossings. We then numerically integrated the one- 
dimensional Schrodingcr equation with the selected potential energy curve and for asymptotically vanishing kinetic 
energy from i?„iin = 40 to i?max = 1940 a.u. The scattering length was obtained from the logarithmic derivative of 
the computed wavefunction at i? = i?max- We verified that our results are converged with respect to variations of the 
integration limits Rmin and i?max- 

We have restricted the lowest magnetic-field strength in Fig. [S] to S = 10 G, because below that the potential well 
starts to disappear. At B = 0, the long-range potential is purely attractive, and additional knowledge about the short- 
range part of the potential would be necessary in order to make quantitative predictions. The strongest field we have 
considered is B — 500 G. The left classical turning point at this field strength is already down to about R — 80. Still 
higher field strengths would again maneuver us into a regime where our theory probably becomes unreliable (see the 
discussion at the beginning of this subsection). However, between 10 and 500 Gauss — and within the approximation 
of pure elastic scattering — our scattering lengths may be regarded as accurate to within about 10%. 

There are two important features about the S-dependence. First, the scattering length vanishes at Bq = 173 G. 
The slope at this point is a'{Bo) = — 0.98/G. Second, a diverges at the resonance field strength Brcs = 339 G, which 
signals the appearance of the first bound state in the long-range potential well. In the vicinity of both Bq and Srcs 
the effective interatomic interaction can be switched between attraction and repulsion. Near B^cs, the ability to tune 
a over an extensive range may be particularly advantageous. 

Interestingly, B^cs and Bq are not independent of one another. We can derive an approximate analytic formula for 
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FIG. 5: Energy expectation values of the s-wave |$o) (Eq. 16711 ) and the five basis states that can couple directly to it (Eq. 1691 '). 
|$i) and |$2) are d-waves; 1"1?3), 1'1?4), and j^s) are g-waves. 



the scattering length as a function of magnetic-field strength: 



a{B) = -—a'iBo)Bo 

OTT 



5tt 



1-tan — \/6/6o- — 



(71) 



which depends only on the parameters Bq and a'{Bo)- Equation (|71|l is based on a semiclassical analysis by Gribakin 
and Flambaum [bJ]. We have employed their analytic result for a hard-core plus a l/_R"-potential. In our case, the 
attraction is basically due to pure dispersion interaction and the repulsion is so strong that replacing it by a hard core 
is not such a bad approximation. In accordance with the discussion following Eq. I|7U|I . we assumed the onset of the 
hard core to scale as B^^^^. The first resonance in Eq. (|71() appears when the argument of the tan-function equals 
7r/2, or B/Bq = (7/5)^ = 1.96. This is in excellent agreement with our numerical data: B^cs/Bq = 339/173 = 1.96. 
The scattering length as predicted by Eq. H71|l . using the parameters Bq = 173 G and a'(Ko) = — 0.98/G estimated 
from our numerical data, is also plotted in Fig. |S1 The analytic formula reproduces the magnetic-field dependence of 
the scattering length remarkably well. 



V. CONCLUSION 



Potential wells associated with ordinary molecular bonds are typically a few electron volts, or 10'' Kelvin, deep; 
the minimum of the potential energy is found at interatomic distances of the order of one Bohr radius. Electronic 
interactions at such short distances are so strong that the energy separation between the ground and excited electronic 
states dissociating to the same atomic configuration is very large in comparison to rotational energies. In a simplified 
picture we might envisage the strong chemical bond as forcing the atoms into a fixed orientation with respect to the 
interatomic axis. Molecular rotation leaves this rigid arrangement virtually unaffected. 

In fact, however, the long-range physics of metastable alkaline-earth atoms is fundamentally different from that 
simplistic picture. If the atoms are separated by a sufficiently large distance, they have constant angular momenta. 
As the atoms approach one another, they experience anisotropic forces that modify the molecular rotational motion 
as well as the relative orientation of the nonspherical atoms. At interatomic distances of a few hundred Bohr radii, the 
energy scales of interatomic interactions and molecular rotations are comparable, and the coupling between atomic 
and molecular angular momentum turns out to be rather efficient. 

The tensorial analysis presented in this paper puts the coupling mechanism into a particularly clear and useful form. 
Our formulation makes maximum use of symmetries, which allows a compact matrix representation of the long-range 
Hamiltonian. It will enable systematic multichannel scattering calculations needed to investigate the role played by 
inelastic collision processes. 

In this study we have numerically diagonalized the Hamiltonian matrix for two metastable strontium atoms and 
obtained adiabatic potential energy curves, the adiabaticity referring only to the distance coordinate R. The curves 
are in general very complicated. We have therefore focussed on one curve — associated with an s-wave channel and with 
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FIG. 6: Coupling strengths of the basis states |<l>i) through j<l>5) (Eq. 1691 1 to the s-wave |$o) (Eq- I67II ). plotted on a logarithmic 
scale. The coupling increases rapidly as the interatomic distance is reduced. It is strongest for the g-waves |$3), |$4), and |$5), 
which can couple to |"I>o) via electric quadrupole-quadrupole interaction. The coupling of |"l?o) to the d-waves |$i) and |$2) is 
mediated by the comparatively weak dispersion interaction. 



the atoms in an experimentally attractive low-field seeking state — that exhibits a novel type of long-range potential 
well. This well arises from a fascinating interplay between the pure l/i?^-attraction in the s-wave channel and the 
strong quadrupole-quadrupole coupling to rotational g-waves attached to lower lying Zeeman manifolds. Recently, 
Avdeenkov and Bohn discovered that the anisotropic forces between polar OH molecules in an electrostatic field lead 
to a similar phenomenon |62|. 

The scattering length in the long-range potential well in metastable strontium can be tuned by varying the external 
magnetic-field strength. We find a resonance at 339 G, which differs quantitatively from the prediction by Derevianko 
et al. |2^], who estimate the resonance location to lie near 1000 G. We believe that the present calculations should 
improve upon the accuracy of the more approximate treatment of Ref. fS^ . 

Finally, the nature of the resonance deserves attention. It indicates the emergence of the first vibrational bound 
state as the long-range potential well under goe s a controlled deformation. This process must be contrasted with the 
physics underlying Feshbach resonances [63, 1641 l65L ISq , where the scattering wave in the entrance channel is brought 
into resonance with a vibrational bound state in a potential well correlating to an energetically higher channel. The 
potential energy curves themselves do not change their qualitative appearance. The strong anisotropic interactions 
in cold gases of metastable alkaline-earth atoms open an entirely new route to shaping potential energy curves and 
creating temporary molecules with bond lengths of the order of a hundred Bohr radii. 
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FIG. 8: Scattering length, a, in the long-range potential well calculated as a function of the magnetic-field strength, B. The 
dots represent numerical data obtained by direct integration of the one-dimensional Schrodinger equation. The solid line is 
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